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ABSTRACT 

The cooler than expected optical-UV transient PSl-lOjh detected by the 
Pan-STARRSl survey is probably related to a tidal disruption event in which 
a He-rich stellar core remnant is implied. The evolution of bound debris during 
the disk phase is studied by solving the hydrodynamic equations. The model 
provides a good fit either of the raising part of the light curve in the bands 
<?pi, t pi s and ip\ or in the early decay. The parameters characterizing this 
optimized model are the mass of the central black hole, i.e., 6.3 x 10 6 M Q 
and the critical Reynolds number 1Z = 10 4 that fixes the viscosity and the 
accretion timescale. Such a high value of 1Z explains the low disk temperature 
and the consequent absence of X-ray emission. The predicted bolometric peak 
luminosity is about 10 45 ergs -1 and the predicted total radiated energy is 
about E rad = 2.67 x 10 51 erg. 
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1 INTRODUCTION 

Different flare events have been detected either in X-rays 
(Bade et al. 1996; Komossa & Bade 1999; Donley et al. 
2002; Esquej et al. 2007, 2008; Cappelluti et al. 2009; 
Maksym et al. 2010) or in the UV region (Gezari et al. 
2006, 2009). These sudden variations in the electromag- 
netic emission of "quiet" galaxies have been interpreted 
as being the consequence of the tidal disruption of a 
star that has passed close to a supermassive black hole 
(SMBH) living in a "dormant" state in the center of 
its host galaxy. For a near parabolic orbit in which the 
periapse distance R p is shorter than the tidal radius 
Rt, tidal disruption is likely to occur and two distinct 
phases in the dynamical evolution of the debris can be 
distinguished. In the first, the "fallback" phase, nearly 
half of the debris is in gas streams having approximately 
ballistic trajectories, which converge and collide at pe- 
riapse, dissipating energy through shocks. Typically, for 
a black hole of mass Mbh ~ 10 7 Mq and a disrupted 
star of m* ~ IM0 the released energy in this phase 
is about 10 erg. Since the dissipated energy in these 
shocks scales as E s h oc M^ 3 m* , in our case, as we 
shall see below, this is expected to be one order of mag- 
nitude less than the "disk-phase. As the bound material 
converges at periapse and dissipates energy, a small ac- 



cretion disk is formed and the second phase initiates as 
soon as the black hole begins to accrete mass, which co- 
incides with the instant of maximum in the light curve. 
This is dictated by the viscous mechanism that controls 
the angular momentum transfer, fixing the timescale in 
which the material is conveyed from the outer regions 
of the disk to the last stable orbit. In this "disk"-phase, 
most of the flare energy is radiated. 

Non-steady accretion disk models have recently 
been built by Montesinos & de Freitas Pacheco (2011a, 
b, hereafter MPUA and MP11B) by solving numerically 
the hydrodynamic and the energy transfer equations. 
These models are able to explain the sudden variation of 
the X-ray emission observed in non-active galaxies that 
are supposed to be related to tidal events like those ob- 
served in NGC 3599 or IC 3599. More recently, Gezari 
et al. (2012) reported an UV-optical flare that has oc- 
curred in a non-active galaxy at z = 0.1696 and which 
was probably originated in a tidal disruption event. This 
flare event discovered in the Pan-STARRSl survey dis- 
plays some particular characteristics like: a) the lack of 
hydrogen lines in the spectra and the presence of a broad 
HeIIA4686 A emission line suggest that the disrupted 
star lost its outer envelope and can probably be iden- 
tified with the He-rich core of an ancient red giant; b) 
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no detectable X-rays emission by Chandra; c) from the 
analysis of the continuum emission, the estimated disk 
temperature is considerably lower than those attributed 
to disks associated usually with X-ray flare events. 

We report here the results obtained from a series 
of models computed with our code, aiming to reproduce 
the observed light curve in different wavelengths and to 
derive the main parameters characterizing the event. 



2 THE MODEL 

The hydrodynamic equations describing the disk evolu- 
tion can be found in MP11A and MP11B respectively 
but, in a few words, we recall here the basic features. 
The code is based on an Eulerian formalism using a 
finite difference method of second-order, according to 
the Van Leer upwind algorithm on a staggered mesh. 
Since the disks considered here are quite small, we have 
adopted an integration grid of 256 ring sectors instead 
of the original 1024 rings adopted in MP11A. The inner 
radius of the grid coincides with the last stable orbit 
while the external radius is defined by the tidal radius 
Rt = iiR*{Mbh/'m*) 3 > where m, and i?» are respec- 
tively the mass and radius of the disrupted star and 
Mbh is the black hole mass. Notice that generally in 
the literature the parameter /j, is taken to be equal to 
unity. In fact, the work by Luminet & Carter (1986) 
suggests that \i = 2.4, which corresponds to the Roche 
limit for a fluid body in a circular orbit. More recent 
investigations of the fly-by problem for a viscous fluid 
body in a parabolic orbit suggest an effective Roche 
limit (R p — R t ) with /i = 1.69 (Sridhar & Tremaine 
1992; Kosovichev & Novikov 1992). In particular, the 
latter authors considered the case in which the disturb- 
ing body is a massive black hole as in the present case. 
Thus, the value [i = 1.69 will be adopted in our grid of 
models. This implies that in our models the most bound 
material returns to periapse after a timescale t m i n that 
is about 4.8 times higher than the case n=l. As we 
shall see later, this explains why Gezari et al. (2012) 
concluded that the disruption occurred only ~ 76 days 
before maximum light while our preferred model indi- 
cates a higher timescale, namely, 1.24 years. The origi- 
nal MP11A boundary conditions were also modified to 
allow an inflow of matter at the external ring, corre- 
sponding to the material that circularizes after shocks 
between converging streams. 

Also as in MP11A, it was assumed that the gravita- 
tional field of the black hole is given by the approximate 
potential of Paczynski-Wiita (Paczynski & Wiita 1980) 
that gives correctly the position of the last stable orbit 
for the Schwarzschild geometry. 

It should be emphasized that in the present model 
the angular momentum transfer is not described by the 
so-called "a" -model introduced by Shakura & Sunyaev 
(1973). The viscosity coefficient is that given by the 
approach of de Freitas Pacheco & Steiner (1976), i.e., 
r\ = 2nrV$/lZ, where r is the radial distance to center 
of the disk, is the azimuthal velocity of the debris 
at that distance and 1Z is the critical Reynolds number 
characterizing the flow. 



Another important aspect concerns the fallback 
rate. According to early investigations (Rees 1988; 
Evans & Kochanek 1989), in a first approximation, the 
fallback rate Rfb is given by the relation 

1^(^) 6/3 , (1) 

where 

t - -2L 5e ( 2 ) 

However, more detailed studies (Lodato et al. 2009) 
based on simulations indicate deviations from eq[T] at 
early phases of the fallback process. In fact, only in the 
late evolutionary stages the fallback rate varies as t~ 5 ^ 3 . 
As in MP11B, we have attempted to include these de- 
viations by modifying eq[T]as 

Kfb ~ 7TX7777 ^ 

[a + [t/tmi n y\ 

In the above equation A is a normalization constant 
(see MP11A for details) and t m m is still given by eq[2] 
Notice that the dimensionless parameter a permits to 
control the instant at which the maximum fallback rate 
occurs. 



3 RESULTS 

A series of models were computed in which the different 
parameters were varied, searching for a best represen- 
tation of data, namely, the light curves of the flare in 
different wavelengths. In all models, we assumed that 
the disrupted star is a He-rich core of 0.23 Mq and a 
radius of 0.33 Rq as in Gezari et al. (2012). Our pre- 
ferred model is characterized by a black hole of mass 
6.3 x 10 6 Mq, a fallback parameter a — 0.06 and a crit- 
ical Reynolds number 1Z — 10 4 . Notice that this value 
is higher than those characterizing our previous mod- 
els (MPUB) but is necessary to get a "colder" disk. 
A higher value of 1Z decreases the local dissipation of 
kinetic energy, producing lower temperatures and con- 
sequently, increasing the energy inflow by advection. 

Contrary what is generally assumed in the litera- 
ture, the beginning of the accretion process by the black 
hole does not coincide with the beginning of the circu- 
larization process or, in other words, with the begin- 
ning of the formation of the accretion disk. This is sim- 
ply because it takes a finite time for the debris to be 
transported from the outer regions down to the inner 
region (last stable orbit) due to the angular momentum 
transfer mechanism. This is clearly seen in fig.l where 
the evolution of the fallback rate is compared with the 
evolution of the accretion rate by the black hole. The 
very initial peak in the accretion rate occurs (for our 
preferred model) about 44.8 days after the onset of the 
circularization process and coincides also with the max- 
imum in the light curve. The mass of the gravitationally 
bound debris is about 0.115 M©, corresponding to a half 
of the mass of the disrupted star, but a small fraction 
(about 15%) is not accreted by the black hole, escaping 
from the system and taking away the angular momen- 
tum stored in the disk. 
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Figure 1 . Evolution of the fallback rate (dashed curve) com- 
pared with the evolution of accretion rate by the black hole 
(solid curve). 



Figure 3. Snapshots of the temperature profile of the disk. 
Labels indicate the corresponding instants of time after the 
beginning of the circularization process. 




Figure 2. Evolution of the average effective temperature of 
the disk. 



Figure 4. Evolution of the bolomctric luminosity. 



The continuous emission of the flare is generally 
fitted by a black-body distribution characterized by an 
effective temperature. In reality, the local effective tem- 
perature varies along the disk surface with a radial pro- 
file that evolves in time. MP11B defined a suitable mean 
effective temperature, which corresponds to an effective 
energy flux that multiplied by the total disk surface 
gives the luminosity. Such a mean effective temperature 
is computed from the relation 



< T 4 >= 



r*P rT 4 (r)dr 

J -H-lso 



(4) 



The evolution of the mean effective temperature is 
shown in fig. 2. The maximum value (~ 39400 K) co- 
incides with maximum light (and maximum accretion 
rate by the black hole) . After approximately 25 days the 
mean effective temperature of the disk drops to about 
25000 K and then decreases more or less steadly at a rate 
of ~ 60K/day. The model predicts at maximum light 



and near the central region a "true" temperature around 
65000 K. One month after maximum, in the inner re- 
gion, the temperature has decreased to about 31000 K 
while in the outskirts of the disk it remains at a value 
around 20000 K. Figure 3 shows different snapshots of 
the temperature profile before, at and after maximum 
light. 

In figure 4, the bolometric light curve is shown. As 
already mentioned, before maximum light the black hole 
is not yet accreting matter and the radiation comes 
from viscous dissipation in the disk that is still being 
formed. Integrating the bolometric luminosity over time 
allows an estimation for the total radiated energy, i.e., 
Erad = 2.67 x 10 erg, in agreement with the lower 
limit estimated by Gezari et al. (2012), based on a black- 
body spectrum characterized by a temperature Tbb > 
30 000 K. Taking into account the total amount of mass 
accreted by the black hole, the mean efficiency for the 
energy conversion is r\ — E ra d/(Mc 2 ) ~ 0.015. 

The resulting spectrum at maximum is shown in fig. 
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Figure 5. Comparison between the "real" spectrum at the 
peak instant with the equivalent blackbody spectra charac- 
terized by the mean temperature of the disk at the peak. 



5 in comparison with the expected spectrum of a black 
body characterized by the mean effective temperature. 
As MP11B have shown, the black body approximation 
underestimates the emitted flux at short wavelengths 
(A < fOOO A) but predicts correctly fluxes at longer 
wavelengths. 

The light curve in filters gpi, rpi and ip\ are shown 
respectively in the upper, middle and lower panels of 
figure 6. Data points are those given by Gezari et al. 
(20f2). According to these authors, who have adopted 
the models by Lodato et al. (2009), systematic differ- 
ences are present between data and theoretical light 
curves during the early rise and the late decay. The 
present model explains quite well both the raising part 
of the light curve and the early decay phase. The former 
is produced while the disk is still being formed and the 
black hole is not yet accreting mass. The early decay, 
just after maximum does not vary as t _5//3 as usually 
claimed but varies faster, decreasing by about one mag- 
nitude per week. The late evolution of the theoretical 
light curve displays the aforementioned behaviour, i.e., 
L oc t~ 5//3 although data in the ip\ filter (lower panel) ~ 
40 days after maximum are slightly below the predicted 
values but still consistent taking into account the obser- 
vational errors. 



4 SUMMARY 

The transient event discovered on 20f0 May 28 by the 
Pan-STARRSf (PSf ) Medium Deep Survey was proba- 
bly originated in a tidal disruption event. According to 
Gezari et al. (20f 2) the disrupted star was a helium-rich 
core resulting from the evolution of a red giant. We have 
used our hydrodynamic code (MPf f A, MPf f B) to com- 
pute the radiation of the bound debris during the phase 
in which a small accretion disk is formed and during 
which most of the flare energy is produced. 

Our preferred model requires that the central black 
hole has a mass of 6.3 x 10 6 Mq almost a factor 2.3 
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Figure 6. Simulated light curves (solid continuous curves) 
compared with data from Gezari et al. (2012). Upper panel 
correspond to gp\ filter, middle panel to rp%, and the lower 
one to ipi . 



higher than that estimated by Gezari et al. (2012) but 
in good agreement with the relation between the black 
hole mass and the stellar mass of the host galaxy. 

No significant X-ray emission related to the flare 
was detected, consistent with the temperature of ~ 35 
000 K estimated by Gezari et al. (2012) that is necessary 
to produce enough UV photons to ionize helium. Our 
model predicts at maximum and near the inner region 
of the disk a temperature of 65 000 K, which is a con- 
sequence of a high critical Reynolds number (TZ = 10 4 ) 
characterizing the turbulent state of the gas. This is 
about 20 times the TZ values defining our previous mod- 
els (MP11B), in which an important X-ray emission is 
present. Notice that in our model the heating rate per 
unit volume e due to viscosity is given by the equation 



27T 2 n 3 

TZ 1 



{ dlgttV 



(5) 



that shows why a higher critical Reynolds number or 
a lower viscosity decreases the heating rate and, conse- 
quently, the disk temperature. 
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The most important difference with the analysis Rees, M.J., 1988, Nature 333, 523 

by Gezari et al. (2012) concerns the chronology of the Shakura, N.I. & Sunyaev, R.A., 1973, A&A 24, 337 

event, which results from the adopted value for the di- Sridhar, S. & Tremaine, S., 1992, Icarus 95, 86 

mensionless parameter /j, defining the tidal radius. If /x ~ 
1.69 and not one, then the disruption of the star oc- 
curred about 1.24 yr before maximum light and the be- 
ginning of the circularization process, about 44.7 days 
before the peak. We emphasize once more that the be- 
ginning of circularization doesn't coincide with the be- 
ginning of the accretion process by the black hole. 

Contrary to the models adopted by Gezari et al. 
(2012), the present model provides a good representa- 
tion of the initial raise of the light curve, when the black 
hole is not yet accreting and the disk is being filling up, 
as well as of the initial light decay, which evolves faster 
than d\og L/d\ogt = 0.6. 

Detailed analyses of tidal disruption events provide 
an excelent tool to study dormant black holes and the 
physics of non-steady accretion disk, in particular of 
the mechanisms reponsible for the angular momentum 
transfer and energy dissipation. 
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